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I. INTRODUCTION 



The commissioning of polarized proton beams with energies from 100 to 250 GeV at 
the Relativistic Heavy Ion Colhder (RHIC) [1] has started a new stage in physics of spin- 
dependent colUsions. In particular, measurement of spin-dependent parton distribution func- 
tions (PDFs) will be the top priority in the experiments with longitudinal polarization. So 
far, these spin-dependent counterparts of the unpolarized PDFs were constrained only by 
the lepton-nucleon deep inelastic scattering (DIS) data [2-10], which concentrates at rel- 
atively small momentum transfers {Q^ ~ 1 GeV^). In contrast, RHIC will explore the 
spin-dependent PDFs in a much larger range of and using a variety of spin-dependent 
particle reactions, including production of Drell-Yan lepton pairs and resonant production 
of VF± and Z° bosons [11]. 

Vector boson production (VBP) with polarized hadron beams, which proceeds through 
annihilation of a quark and antiquark at the Born level, is the most natural candidate to 
probe spin-dependent quark polarizations [12-14]. Furthermore, in pp collisions the Born 
level cross sections are sensitive to the distributions of sea quarks. In this sense, VBP 
complements production of jets, pions, and heavy quark flavors, which primarily probe the 
gluon distribution. No wonder that extensive bibliography is dedicated to spin-dependent 
production of Drell-Yan pairs [15-28] and massive electroweak bosons [26, 29-36]. 

The production of bosons presents a particularly interesting opportunity to learn 
about the quark spin structure due to the maximal violation of space-reflection parity in the 
qqW coupling and non-trivial mixing of the quark flavors through the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix [37, 38]. While the former feature allows non- vanishing single- 
spin cross section asymmetries (which are simpler than the parity-conserving double-spin 
asymmetries), the latter feature facilitates the study of the flavor dependence of the sea 
quark PDFs. The issue of flavor symmetry breaking in the polarized quark sea was brought 
in the limelight recently by the results on semi-inclusive hadroproduction in spin-dependent 
DIS by HERMES Collaboration [39]. In that measurement, the issues of higher-order and 
higher- twist corrections are an important consideration due to small values of Q^. At the 
same time, W boson production at RHIC has a potential to pin down the quark sea flavor 
structure at large Q^, i.e., at energies where perturbative quantum chromodynamics (PQCD) 
is truly valid. 
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The original method [11, 40] for the measurement of polarized PDFs in W boson pro- 
duction at RHIC relies on the reconstruction of the spin-dependent distribution dcr/dy over 
the rapidity y of the W bosons. The method is based on the observation that the Born 
level single-spin asymmetry of da/dy is described by a simple theoretical expression, which 
further reduces to the ratio Aq{x)/q{x) of the polarized and unpolarized PDFs when the 
absolute value of y is large. Unfortunately, the application of this method is obstructed 
by specifics of the detection of the W"^ bosons at RHIC. Since neither of RHIC detectors 
controls the energy balance in particle reactions, the energy and momentum of the neutrino 
in the decay are completely unknown, so that the information about the momentum of the 
W boson is incomplete. As a result, the determination of the rapidity y of the W boson 
from the observation of just one charged lepton is generally impossible. 

It can be shown that the ambiguity in the reconstruction of y reduces to the uncertainty 
in the knowledge of the transverse momentum qx of the W boson [11]. If qx were known 
exactly, the rapidity y of the W boson can be derived (up to a two-fold ambiguity) from 
the measured rapidity ye and transverse momentum pxe of the charged lepton (see the 
accompanying paper [41] for more details on this reconstruction method). For instance, if 
all W bosons were produced through the Born process qq ^ W, the transverse momentum qx 
would be zero, and the correct solution for y can be chosen statistically for the leptons with 
large absolute values of the lepton rapidity y^, i.e., escaping near the beam pipe direction. 

In reality, the bosons carry a non-zero transverse momentum due to QCD radiation, 
with the most probable magnitude in the unpolarized case of about 2 GeV. Most of the 
bosons obtain a non-zero qr through the radiation of soft and collinear partons, which 
cannot be approximated by finite-order perturbative calculations. Since the power series in 
the strong coupling as does not converge in the small-g-r region, summation of dominant 
logarithmic terms through all orders of this series is needed. This all-order summation can 
be realized with the help of the methods that were studied in a substantial detail [42- 
58] in the unpolarized case. Note, however, that the properties of the multiple parton 
radiation depend on the spin of the initial hadrons. The spin dependence of the collinear 
radiation off the external parton lines can be seen from an explicit calculation. An additional 
spin dependence can be contributed by the unknown nonperturbative dynamics of strong 
interactions at large distances. Thus, conclusions about the nature of the multiple parton 
radiation in the polarized hadronic collisions cannot be inferred from the unpolarized case. 
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and an additional study is required to estimate the spin dependence of such radiation under 
the RHIC conditions. 

The main goal of this paper is to provide a theoretical framework for such a study. We 
present a complete formalism for resumniation in VBP with the proton beams of an arbitrary 
longitudinal polarization. Furthermore, we explicitly account for the decay of the vector 
bosons into the observed leptons, which is important because of the complicated geometry of 
RHIC detectors. The results are derived at a one-loop level of PQCD, and the summation of 
large logarithms is performed with the help of the impact parameter resummation technique 
[47]. The detailed discussion of this lepton-level resummation formalism in the unpolarized 
case can be found in Ref. [50]. The theoretical results are presented for arbitrary couplings 
of the vector boson, which makes the results of this paper also applicable to lepton pair 
production mediated by photons or bosons. Our eventual product is a numerical program 
for Monte-Carlo integration of fully differential resummed cross sections in the presence 
of experimental cuts. In the accompanying paper [41], we present numerical results for 
the single-spin asymmetries in W boson production and estimate the sensitivity of these 
asymmetries to different models of polarized PDFs. 

To review the already available literature, we note that the finite-order asymmetries of 
order 0{as) in the Drell-Yan process are currently available for distributions in full lepton 
pair momentum [19], invariant mass [22], and rapidity [23]. Analogous distributions for W 
boson production were published in Refs. [30, 35, 36]. Recently, the fully differential 

distributions of the Drell-Yan pairs were obtained in Ref. [27]. Furthermore, Refs. [19, 30] 
presented the 0{as) resummed single- and double-spin cross sections in the narrow width 
approximation, and the 0{as) Sudakov factor was explicitly demonstrated to be independent 
from the polarization of the proton beams. Note that most of the above publications (with 
the exception of the forward-backward asymmetry in Ref. [35] and polar angle distribution 
of the Drell-Yan leptons in the lab frame in Ref. [28]) discuss the behavior of the whole 
lepton pair, rather than the distributions of the individual leptons. 

In this paper, we go beyond the previous results in several aspects. First, we present the 
fully differential finite-order and resummed cross sections for the decaying electroweak vector 
bosons, i.e., we completely account for the spin correlations between the hadronic subsystem 
and leptonic final state. In other words, the cross sections in this work are fully differential 
in the individual lepton momenta, rather than in the momentum of the lepton pair. These 
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cross sections are implemented in a Monte-Carlo integration program. Such lepton-level 
analysis requires calculation of several additional angular structure functions, which do not 
contribute in the narrow width approximation. Furthermore, the Qt resummation is needed 
not only for the parity-conserving angular function 1+cos^ 6', which contributes to the boson- 
level cross section, but also for the parity- violating angular function 2 cos6', which affects the 
angular distributions of the decay products. The calculation for the parity-violating term 
2 cos 6 is substantially more complex, because it involves a large number of 75 matrices (and 
Levi-Civita tensors) coming from both spin projection operators and axial couplings of the 
electroweak Lagrangian. In contrast, the parity-conserving term 1 + cos^ 6 depends on 75 
matrices and Levi-Civita tensors from the spin projection operators only. 

It is well known that special care is needed to deal with the 75 matrices in n 7^ 4 
dimensions. In order to tackle this issue efficiently, we have evaluated the n-dimensional cross 
sections using the dimensional reduction method [59-61]. We then converted our results to 
the conventional MS factorization scheme, which utilizes dimensional regularization [62, 63]. 
From our calculation, we have found that the spin-dependent resummed cross-sections in 
[19, 30] could not be used with the existing parametrizations of the polarized PDFs in the MS 
scheme, due to the different factorization scheme used in those papers. In addition, Ref. [30] 
has used an unconventional normalization for the n-dimensional single-spin cross section in 
the qG channel. The present work has paid a special attention to rectify those inconsistencies 
and obtain the hard cross sections compatible with the existing phenomenological PDF 
parametrizations in the MS scheme. 

The paper has the following structure. In Section II, we introduce the notations for the 
kinematical variables and spin-dependent cross sections. Section III discusses the regular- 
ization of cross section singularities by continuation of observables to n 7^ 4 dimensions, as 
well as the treatment of the 75 matrices. In Section IV, we present in detail the 0{as) 
finite-order and resummed cross sections. Section V contains the summary of our results 
and conclusions. The appendix contains the expressions for all structure functions that 
contribute to the lepton-level 0{as) cross sections. 
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Figure 1: Notations for the momenta and helicities of the participating particles. 



II. NOTATIONS 



A. Kinematical variables 

We discuss inclusive production of a lepton pair A{pa) + B{pb) ^*('?) + X ^ h{h) + 
^2(^2) + X mediated by a virtual vector boson V. Our notations for the particle momenta 
and helicities in this process are shown in Fig. 1. 

The momentum g^' of V can be conveniently parametrized by the invariant mass Q, 
rapidity y, and transverse momentum qt of the vector boson in the laboratory frame. We 
have 

Q' = i''q„ (1) 

Qt = -QtQt,,, (3) 

where is the component of the momentum g'^ that is orthogonal (in the relativistic sense) 
to the momenta and p'^: 

- - y^^PS - /^rt. (4) 

[Pa ■ Pb) [Pa ■ Pb) 

For any vector n^, the "plus" and "minus" components are defined as n± = (n° 
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It will be useful to introduce the hadronic and partonic Mandelstam invariants, defined 

by 

s = {PA + Psf; (5) 

t = {PA- qf\ (6) 
u = {pb- q?; (7) 

S = {Pa + Pbf] (8) 

t= {Pa-qf\ (9) 
u = {p,- qf. (10) 

A caret denotes quantities at the parton level. They are defined in terms of the parton 
momenta = ^aP^Ai Ph = ^bPsi where and are the light-cone momentum fractions, 
which satisfy the constraints x^i < < 1, < < 1, with xa,b = Qe^^ l\fs. Throughout 
the discussion, all hadron masses are neglected. 

Our discussion will use two coordinate frames. The first frame is the laboratory frame, or 
the center-of-mass frame of the initial proton beams. The second frame is a special rest frame 
of the produced vector boson (CoUins-Soper frame [64]), in which the axis Oz bisects the 
angle between the vectors p^^ and — ps, where Pa are pb are the initial hadrons' momenta 
in this frame. The coordinate transformation from the lab frame to the Collins-Soper frame 
involves a boost in the direction of motion of V and a rotation around the Oy axis (see 
Appendix A in [50] for the explicit coordinate transformation matrix). The components of 
the momenta p\ and in the Collins-Soper frame are 



where Mr is the transverse mass of the vector boson, Mt = a/Q^ + q^- The momenta 
and ^2 of the final-state leptons in this frame are 

V{ = Y sin 6* cos sin 6* sin y?, cos6'^ , (12) 

I2 = ^ ( 1, — sin 6' cos V9, — sin 6' sin — cos6']. (13) 
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Instead of using and I2 directly, it is convenient to operate with their linear combina- 
tions, 

q'' = l',+lt^=(Q,0,0,0] (14) 



and 

1^2 = li — I2 = Q ^0? sin 6 cos sin 6 sin y9, cos 6^ . (15) 

While describes the motion of the lepton pair as a whole, 1^2 specifies the motion of the 
individual leptons in the decay of the vector boson. To separate the dynamics of vector 
boson production from the dynamics of vector boson decay, we decompose the lepton-level 
cross section in a sum over the functions Ap{9, (p) of the angles 6 and (f in the Collins-Soper 
frame: 

lirik^-t'nQ\y.4)A,(e.^). (le) 

In this equation, dQ = d cos 9 d(p denotes an element of the solid angle in the Collins-Soper 
frame. ''T{Q'^, y, q^) is the structure function corresponding to the angular function Ap{6, ip). 
For an arbitrary chirality of the electroweak couplings, VBP at 0{as) receives contributions 
from six angular functions Ap{9,ip)-} 

A^i = I + cos^ 6; Ao = ^{l-3cos^9); 

Ai = sin 29 cos ip; A2 = - sin^ 9 cos 2ip; 

As = 2 cos 9; = sin6'cosv3. (17) 

The angular function A_i{9,(p) is also sometimes denoted as Lq{9,(p) (see, for instance, 
Ref. [50]). It will be shown below that the structure functions ~^T, °T, ^T, and are 
generated by the vector part of the electroweak current. On the other hand, the structure 
functions and are generated by the axial part of the electroweak current. This feature 
leads to the following important consequence. When the cross section is integrated over the 
complete solid angle Q of the lepton decay (as, e.g., in the calculation of a boson-level cross 
section), all angular functions except A_i are integrated out. The resulting cross section 
is sensitive only to the vector part of the electroweak current, so that, for instance, the W 
boson cross section can be straightforwardly derived from the Drell-Yan pair cross section. 
On the other hand, if one is interested in the angular distributions of the final-state leptons, 



^ At 0{ag), there will be additional contributions to the longitudinally polarized cross section, which are 
proportional to sin 20 sin 1^9 and sm^29sm(f [65]. These contributions lead to non- vanishing single-spin 
asymmetries at 7^ in the lepton distributions in the parity-conserving case [65-67]. Such contributions, 
however, do not appear at the order 0{as) discussed here. 
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the axial part of the electroweak current cannot be ignored. In that case, additional care 
is needed in the calculation of the parity-violating structure functions and ^T, which 
are affected by the 75 matrices from both the spin projection operators and electroweak 
couplings. 



B. Polarized cross sections 



We now introduce the special notations that will allow us to write the cross sections for 
arbitrary proton polarizations in a compact form. Let the polarized protons A and B have 
helicities Ha and Hb, respectively. The allowed values are +1 for the right-handed helicities 
and —1 for the left-handed helicities. The cross section for this combination of the helicities 
is denoted as a^^^B _ Yoy brevity, the superscripts of a^^^'-B will only show the signs of the 
helicities; that is, cr+^'+^ = cr++, a^^'^^ = (t+'^b ^ etc. 

The following combinations of the helicity cross sections will be called "the cross sections 
that are unpolarized (U) or polarized (P) on the side of the proton A": 

Similarly, "the cross sections that are unpolarized (polarized) on the side of the of the proton 

5" are 

^{h^m ^ ^h^+^^h^-^ (20) 

and 

^ihA,P) = ^hA+ _ ^hA-^ (21) 

respectively. 

Next, we introduce the unpolarized, single-spin, and double-spin cross sections a^'^''^\ 
and a^^'^^: 



= \ E (22) 

hA,hB=±i 

= i E hAa'-'^-; (23) 

= J Yl hAhBa'^^^^. (24) 
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hA,hB=±l 



The normalization factor 1/4 in (22-24) corresponds to the number of the allowed helicity 
combinations for the initial-state massless hadrons. The single-spin cross-section a^^'^^ 
in Eq. (23) corresponds to the polarized beam A and unpolarized beam B. The fourth 
independent linear combination of the helicity cross-sections, the single-spin cross-section 
a^^'^\ can be obtained from the single-spin cross-section a^^'^^ by interchanging the indices 
of the proton beams, A ^ B. At 0{as), the single-spin cross section (23) is non-zero only 
if VBP violates parity with respect to the spatial reflection. In other popular (but less 
uniform) notations, a^'^'^\ a^^'^\ and a^^'^^ are denoted as a, ALa^ and I^llo, respectively. 
Hence, the conventional single- and double-spin asymmetries are constructed as 

= (26) 

where the deflnition for the single-spin asymmetry Al does not include an additional minus 
sign. This choice leads to the negative values of Ai for the left-handed interactions of 
fermions at the Born level, and it is the same as the deflnition in Refs. [35, 36]. 

The QCD factorization expresses the hadron-level cross sections (J^^'^ in terms of the 
"hard" parton-level cross sections d^ab'^''{fJ'F) and parton distribution functions: 

^A^'"= E E /'^^<^/'^e^^a6'H/^F)A.//..(ea,/iF)/w/^,(e6,/^F). (27) 

a,b=G, u , d 

Here fha/hAi^a, fJ'p) denotes a helicity-dependent parton distribution function (PDF), i.e., 
the probability of flnding a parton a with the momentum = ^aPA and helicity ha in a 
hadron A with the momentum and helicity Ha- The parton-level cross section o^ab^^'ilJ'F) is 
separated from the PDFs fha/hAi^a, fJ^p) at a factorization scale fip, which in our calculation 
is assumed to coincide with the QCD renormalization scale. 

Since the strong interactions conserve parity, only two linear combinations of the helicity- 
dependent PDFs fha/hA^ia^ I^f) are independent. These combinations will be called the 
unpolarized and polarized PDFs, respectively: 

ff/li^a^ f^p) = f+/+i^a, fip) + f-/+{^a, f^p) = f+/-{^a, fip) + f-/-i^a, f^p), (28) 
fi^li^a, f^p) = f+/+{U, f^p) - f-/+{^a, f^p) = -/+/-(^a, f^p) + f~/~{U, f^p)- (29) 

In common notations, fa^l{^a,fJ'p) = fa/AiC.a, fJ'p), and fa^l{C.a, f^p) = Afa/A{C.a, f^p)- When 
written in terms of the unpolarized and polarized cross sections, the factorization formula 
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Table I: The couplings of the electroweak vector bosons V to the fermions. The left couplings 
are and the right couplings are fR,gR] see Eqs. (31) and (32) for their definition, g is the 

electroweak coupling, aEMifJ-) is the running fine structure constant, and — sin^-j^ — cos^^) 
is the sine (cosine) of the weak mixing angle. Qf is the fractional charge of the fermions [Qu = 2/3, 
Qd = —1/3, Qu = 0, Qf,- = — 1). T3 is the eigenvalue of the weak isospin for the fermion 
[T^ = T- = 1/2, Ti = Tf =-1/2). 

V fL,gL fn, gR 

7 ^y/AnaEMilJ-) ^Y^47ra£;A/(/x) 

g/{2V2) 
^° g{n-Qfsl)/{2c^) -gQfsl/{2c^) 



(27) sums only over the types of the partons and not over their helicities: 



'^Tb'''''^ = 1^ I dial dib^T''''\f^F)f:;f{u,f^F)n'/]!>{ib,f^F), m 

a,b=G, u , d 

v\fhere Ha, Hb denote unpolarized or polarized quantities: 

Ha,b = U, P. 

Eq. (30) already illustrates the usefulness of the indices U and P: in one equation, it covers 
the factorized representations for all three cases of the unpolarized, single-spin, and double- 
spin cross sections. In the following parts of the paper, many expressions will be presented 
in this compact and uniform notation. 



C. Constant overall factors 

The results in this paper cover production of virtual photons, W, and Z bosons. We 
use the on-shell scheme and improved Born approximation to parametrize the electroweak 
parameters in the considered processes. Let fi, Jr and gi, Qr be the chiral couplings entering 
the leptonic vertex ^1^2^ and quark vertex qjQkV, respectively. Here the quark flavor indices 
are j = m, d, s, ... and k = u,d, s, ... .The corresponding interaction Lagrangians are expressed 
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as 



Ai/.y = - 75) + + 75) J , (31) 

and 

^^imv = [adl - 75) + gnil + 75)) Vfk, (32) 
respectively. The matrix Vj% defines the flavor structure of the qj^V vertex and is given by 



, CKM matrix elements for W^, 

yfk = y-k, = \ ^ _ . • (33) 



6,k for 7%Z0. 



The values of the leptonic chiral couplings /l, fji and quark chiral couplings ql, to 7*, 
^ and bosons are specified in Table I. They are expressed in terms of the running fine 
structure constant aEAiifJ'), weak coupling g, and sine of the weak angle sinOw, calculated 
as 

/ = 4V2M^Gf (34) 
and 

sin^Ow = 1- M^/Ml (35) 

from the input values of the Fermi constant Gp, W boson mass My/j and boson mass 
Mz- The nice feature about this parametrization of the electroweak couplings is that the 
higher-order electroweak radiative corrections in VBP are reduced. 

Various constant factors in front of the cross sections will be absorbed in the hadronic 
normalization constants and leptonic normalization constants af: 



'ah 

27r 



(fi'L ifi'fl) iKbP, a,b = G, u, u, d, d,. . . ; (36) 



a 



± 



1 Q 



2 



{fl ± fl) ■ (37) 



" 327r2(Q2_M2)2 + r2.QV^v 
In these equations, = 3 is the number of colors, My and Ty denote the mass and width of 
the vector boson. The indices a and b denote all possible active parton flavors, including the 
gluons. When both a and b correspond to the quarks or antiquarks (a, b = u, u, d, d, . . 
the matrix element Kfe is given by Eq. (33). If one of the indices in Vat corresponds to the 
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a) b) c) d) 




e) f) g) h) 



Figure 2: Partonic subprocesses contributing to vector boson production. 

gluons, we implicitly assume summation over the quarks or antiquarks on the side of the 
gluon: 

\yjG?= i^i^i'' \^ok?= E i^^-^i'' ^tc. (38) 

k=u,d,... j=u,d,... 

Unless explicitly stated otherv^^ise, the subscripts j and k run over the quark and antiquark 
flavors, respectively (j = u,d, and k = u,d, ...). They should be distinguished from the 
indices a and b, which can take the value of any parton flavor (a, b = G,u,u,d,d, ...). 

III. REGULARIZATION OF SINGULARITIES: GENERAL PROCEDURE 

The complete set of Feynman diagrams contributing at order 0{as) is shown in Fig. 2. 
The calculation involves cancellation of ultraviolet singularities in the sum of all virtual 
diagrams (Figs. 2b-2d), cancellation of soft singularities in the sum of virtual and real emis- 
sion corrections (Figs. 2b-2h), and factorization of collinear singularities in real emission 
corrections (Figs. 2e-2h). Those singularities should be flrst exposed by intermediate regu- 
larization, which can be achieved by continuation to n = A ~ 2e dimensions. In our case, 
this approach involves the necessity to define 75 matrices in n 7^ 4 dimensions. As is well 
known, such continuation cannot be mathematically consistent and preserve symmetries of 
the classical Lagrangian at the same time. For instance, dimensional regularization (DREG) 
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[62, 68, 69] assumes that the 75 matrix in n > 4 dimensions is a purely four-dimensional 
object satisfying the following commutation relations: 

{75, 7a.} = 0, if /i < 4; [75, 7^] =0, if /i > 4. 

Since the conventional MS scheme utilizes dimensional regularization, we will sometimes 
refer to this scheme also as the "DREG factorization scheme". The DREG choice is math- 
ematically consistent [63]. At the same time, it requires to treat vector components in the 
(n — 4)-dimensional subspace differently from vector components in the four-dimensional 
subspace. The proliferation of extra (n — 4)-dimensional terms considerably complicates the 
algebra and, more importantly, leads to the violation of chiral symmetry [70] and supersym- 
metry [71, 72]. In spin-dependent QCD, this implies non-conservation of the quark helicity 
in the process of gluon radiation. Those symmetries have to be restored order by order by 
introducing additional counterterms. At order 0{as)i a well-known example of this feature 
is provided by an additional renormalization of the one-loop virtual QCD corrections needed 
to restore the Ward identity for the axial electroweak current [70]. Another relevant exam- 
ple is given below in the discussion of Eq. (57), where an additional finite renormalization 
is performed in DREG in the part associated with the quark-gluon splitting in 7^ 4 di- 
mensions. In both cases, the additional renormalization restores conservation of helicity in 
the processes with radiation of soft or collinear partons, which is otherwise violated by the 
evanescent {n — 4)-dimensional terms appearing in the DREG calculation. 

On the other hand, alternative approaches sacrifice certain features of the full theory, such 
as the cyclic permutability in Dirac traces [73] in the anticommuting 75 scheme [74, 75]; or 
they are mathematically inconsistent at higher orders of the perturbative expansion. In a 
popular alternative to DREG, dimensional reduction (DRED) [59-61], the spinor indices are 
kept in four dimensions, while particle momenta are declared to have n < A components. 
By first evaluating the Dirac traces in four dimensions and then evaluating loop integrals 
in n dimensions, one explicitly preserves the quark helicity in the QCD vertex. Since the 
DRED method treats all spin components on the same footing, it is also algebraically simpler 
than the DREG approach. Nonetheless, contradictions in the DRED framework are present 
in the diagrams with more than two loops [76, 77]. These contradictions, which reflect 
inconsistency in the decomposition of the 4-dimensional momentum space into 4-dimensional 
subspace (associated with the spin indices) and (4 — n)-dimensional subspace (associated 
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with the momentum indices), do not occur at lower orders of the perturbative expansion. 
Hence, it is justified to use the DRED method at the one-loop order of PQCD. 

Due to the large number of 75 matrices to be handled, and to avoid entirely the issue 
of additional symmetry-enforcing renormalizations in the DREG scheme, we have chosen 
to perform the current one-loop calculation in the DRED framework. In particular, soft or 
collinear radiative corrections in DRED automatically preserve the chirality structure of the 
electroweak vertex, so that such corrections do not affect the angular distributions of the 
leptons observed at the Born level (cf. further discussion in section IV). Once the DRED 



results are available, they can be transformed to the MS (DREG) factorization scheme 
to make them compatible with the existing phenomenological PDF parametrizations. The 
detailed description of this transformation can be found in Refs. [22, 35, 60, 61, 78], and 
we refer the reader to those papers for more information. In the following, we will outline 
only those aspects of this method that are relevant to the derivation of the resummed cross 
sections. 

We start by calculating the parton-level cross sections 

,Z^?'' = ^77^5 [s + t + u- Q'] H,^L^\ (39) 

In this equation, the (purely four-dimensional) leptonic tensor L^^ is given by 

L'"" = 4| [fl + fl] i-q'g^'' + q^q"" - f^.l^,) +i[fl- fl] e'^'^'^^ q^{h,) , (40) 

and the hadronic tensor H^y includes squared matrix elements corresponding to the Feynman 
diagrams in Fig. 2. We first evaluate the Dirac traces with the help of the package TRACER 
[79] and keeping all particle spins in four dimensions. The projections on the unpolarized 
(f/) and polarized (P) initial parton states are realized by inserting an appropriate spin 
projection operator 

Vf' = \l. or ) = ^75 (41) 
on each incoming quark leg, and 

if a 13 , P a\ 
l-^G ~ 2 I f I ' ~ f Pc'iPh'S [4:2,) 

on each incoming gluon leg. 
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We then contract the resulting four-dimensional tensors with n-dimensional particle mo- 
menta, neglecting the mismatch between the four- and n-dimensional indices. We find that 
the individual virtual diagrams (Figs. 2b-2d) contain ultraviolet singularities, which, how- 
ever, cancel in the sum of all virtual diagrams. Furthermore, the real emission contributions 
contain the singularities when the unobserved final-state parton is soft or collinear to one 
of the incoming hadrons. In the soft and collinear limits the transverse momentum qr ap- 
proaches zero. The soft singularities (appearing as the poles proportional to l/e^ and 1/e) 
cancel in the sum of all 0{as) diagrams. The remaining task is to factorize the collinear 
singularities according to the factorization scheme that was used to define the hadron-level 
PDFs. 

It can be shown that the choice of the 75 prescription affects only functions C^^^{^,h, jj) 
(where H = U ov P) m the resummed cross section (see the detailed discussion of these 
functions in the next section). The C-functions contain (n — 4)-dimensional parts of the 
n-dimensional splitting functions, which are different in the DRED and DREG schemes. 
An 0{as) n-dimensional splitting function -Pa/^'''"(0 can be calculated in the same way as 
the conventional four-dimensional splitting function -Pj/^ (Oi except that it also retains the 
terms of order 0{e). Hence, it contains an additional finite part eP^/i, (Oi which vanishes 



as n — > 4: 

= i^i/?(0 + ^i^i/?'^(0; (43) 

= i'i/?(0- (44) 
The four-dimensional splitting functions -P^/^HO needed for our calculations are 

Pl'/iiO = i + (1 - 0") («) 

in the unpolarized case [80-82], and 

n'/';'K) = c^-(^)^=CK). («) 

P^iS) = 2^-1 (48) 



in the spin-dependent case [82, 83]. In the unfactorized parton-level cross section (39), the 
functions -P|/^^'"(Oi which are implicitly contained in the hadronic tensor Hf^^, are multiplied 
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by the collinear poles i.e., they appear in the terms 



as 1 piH) 



27re 



27r 



(49) 



The first piece on the right-hand side of Eq. (49) is an 0{as) contribution to the parton-level 
distribution function fl^Ji^, s)- Hence, in accordance with the factorization theorem (30) for 
the case when A = a', and B = b', it is factorized out of the parton-level cross section aa'b' 
to obtain the hard section aa"b"- On the other hand, the second piece in Eq. (49) remains 

(H) £ 

in the hard part. As a result, the terms -P^/fe (0 included in the 0{as/TT) coefficients 
C^^^^^\^,bfi) of the C-functions, as 



C, 



(50) 



and 



C. 



(mm 

i/G 



jjh 



i,j = u, u, d, d, 



(51) 



The meaning of the constants Ci, C2, and 60 is discussed in the next section. 

Since in the DRED scheme the particle spins are exactly four-dimensional, the {n — 4)- 
dimensional parts -Pj/^ (0 of the DRED splitting functions identically vanish: 



DRED: Pi/^)'^(0 = 0. 

Therefore, in the DRED scheme the C-functions (50) and (51) reduce to 



(52) 



DRED: C, 



(e, bfi) 



- CF5{l-i) 



19 _ 
L16 4 
fib 



In' 



-3^ 
boC2 



(53) 
(54) 



In the DREG scheme, the (n — 4)-dimensional parts of the unpolarized [84] and polarized 
[85] splitting functions are given by 



DREG: P^^]''{0 



<?/<? 

q/G 



+ -6(1-0 



(0 = e'-e, 



(55) 
(56) 
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g/G 



(0 



(0 



Cf 
Cf 



3(1-0 + 2^(1 + 



-(1 

1. 



p 



{U),e 



(0, 



(57) 
(58) 



In the quark-initiated polarized splitting function Pg^^'^^iO (see Eq. (57)), an additional 
finite renormalization term Zi;{^) = —4(7^(1 — was included to restore the quark helicity 
conservation in the quark-gluon vertex, in accordance with the existing definition in the 
DREG scheme [86-88]. With the help of Eqs. (50) and (51), the C-functions in the DREG 
scheme are found to be 



DREG: Cg'I'iK.V) = C'™((M = 



23 
16 



- CF6il-0 

C'^'(e,¥) = ^e(i-o-Ooin( 



^ 1 2 

4 



fib 



boC2 



and 



C. 



(p)m 

i/G 



(e, hfi) 



^(i-o-nTc(^)in( 



(59) 
(60) 

(61) 



To summarize this part, the difference between the DRED and DREG factorization 
schemes is contained entirely in the functions Cl^^{^, b, fi) of the resummed cross section. The 
form of the scheme-dependent terms in the C-functions is determined by the n-dimensional 
splitting functions in each factorization scheme. Hence, to transform the results from the 
DRED scheme to the DREG scheme, one simply finds where the DRED (n — 4)-dimensional 

(H) € 

functions P^/b (0 would appear in the C-functions if they were not identically zero. One 

(H) 6 

then inserts at those places the DREG functions P^/i, (0 given in Eqs. (55)- (58). 
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IV. CROSS SECTIONS 



A. Finite-order results 



Let us now explicitly present the VBP cross sections, starting from the lowest order (LO) 
hadron-level cross section (Fig. 2a): 



da 



(Ha,Hb) 
AB 



dQ'^dyd'^qj-dQ 



LO 3=u,d,- k=u,d,. 

+ (j ^ k) 



q{Ha,Hb)i 
jk 



(Hb), 



(62) 



As discussed in Section II, Ha and Hb can take one of the two values {U or P) corresponding 
to the cross section that is unpolarized (or polarized) on the side of the respective initial-state 
hadron. The index j runs over the active quark flavors, while k runs over the active antiquark 
flavors. The functions G^^'^'^^\6,Lp) in Eq. (62) are composed of the normalization factors 



CT; , a.j^ from Eqs. (36), (37) and angular functions A_i{6,ip) and A2,{6,(^) from Eq. (17): 



jk 



]k 



(63) 



The flavor-space tensor e^-^ determines the change of the sign under the interchange of quarks 
and antiquarks: 



1, 



^ki 



-1. 



(64) 



Note that the parity- violating angular function A?, = cos 9 does not contribute if V couples 
to the fermions through the vector current, i.e., if cr^ = crT^ = 0. From the last equation (63), 
we see that the double-spin parton-level cross section is equal to the unpolarized cross section 
with the minus sign. This equality is a direct consequence of the chirality conservation in 
the electroweak couplings (31) and (32). 

Next, consider the 0{as) virtual corrections shown in Figs. 2b-2d. In n = 4 — 26 dimen- 
sions, the contribution of the virtual corrections is given by the LO cross section scaled by 
a function $(05, e): 



da 



AB 



dQ^dyd^'-'^HTdn 



<l>{as,e) 



da 



{Ha,Hb) 
AB 



0{as), virtual 



dQ'^dyd^-'^qrdn 



(65) 



LO 
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as /47r/i2,y 1 / 2 3 2 



where 

The finite part of $(a5,e) depends on the prescription for the 75 matrices in n 7^ 4 dimen- 
sions. In the DRED scheme 6y = 1, while in the DREG scheme 6^ = 0. 

The remaining contributions are from the 0{as) real emission processes shown in Figs. 2e- 
2h. The real emission cross section can be expressed in terms of the PDFs f^^\^a, Pf), 
angular functions Ap(6', (/?), and parton-level structure functions '^Tab{HA, HB,'s,t,u): 



{Ha,Hb) 



2^ / ^ / A ^S{s + t + u-Q 



2tt'^S ^ _ Jq Jo S 

0{ag),real a,b=G,u,u,d,d,... 



dQ'^dy dq^dVL 



In Eq. (67), Glflb '^'^^ the spin-dependent combinations of the constant factors af and 



^ p,ab 



CT^j defined in Eqs. (36) and (37): 



MU,U) ^ MP,P) _ ) ^i^ab for P = -1, , 1 , 2; 

^p,ab ^ p,ab 



^e^ab for P = 3,4; 

(^i^^ab for P = -1,0,1 
^i^lb for P = 3,4. 



MP,u) ) ^I^ab for p = -1,0,1,2; 

^P,ab = S , „ ^ . (68) 



Csu{3) is a color group factor, which is equal to Cp = {N^ — l)/{2Nc) = 4/3 in the annihi- 
lation subprocess qq GV and Tr = 1/2 in Compton scattering subprocess q G ^ q V. 
The structure functions ^TabiHA, HB,'s,t,u) in the unpolarized case {Ha = Hb = U) can 
be found in Ref [50]. Our new result is the calculation of these functions in the polarized 
case {Ha = P and/or Hb = P)- The explicit expressions for these functions are given in 
the appendix. 

For any polarization of the proton beams, the structure functions ~^Tab{HA,HB,^,t,u) 
and ^Tab{HA, HB,'s,t,u) have a singularity at qx = 0, which corresponds to the emission of 
soft or collinear real partons. The leading-logarithmic part of the 0{as) finite-order cross 
section {asymptotic piece) at small, but non-zero (0 < ^ Q) is 



dQ'^dy dq^dVt 



s 7T 27rqi 

O{as),qT-^0 j=u,d,- k=u,d, 
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X 



Cp In 



(1t_ 



{Ha), 
I A ' 



XA,f^F)fi^B^{xB,fJ'F) 



(69) 



Here the functions G^^^'^"^^ {6, if) are the same as in Eq. (63), and Pa^f,^\0 the one-loop 
splitting functions in Eqs. (45)- (48). The convolution of two functions fa/b{x) and gb/c{x) is 
defined by 

[fa/b ® gb/c]{x) ^ I jfa/b{j)gb/c (0 • (70) 

b=G,u,u,d,d,... 

This definition includes summation over the repeating parton index h. 

In accordance with the Kinoshita-Lee-Nauenberg theorem, the soft singularity, which 
appears in the cross section (69) at = 0, is canceled when this cross section is added to 
the 0{as) virtual contribution (65). Furthermore, the collinear singularities are absorbed 
in the parton distributions. The complete small-gy cross section (including both real and 
virtual corrections) is given by 



da 



{Ha,Hb) 
AB 



dQ'^dy dq^dfl 
i5^(q^)(G(.f'^-He,^) 







O{as),qT^0_ 


real+virt. 



da 



AB 



dQ^dy dq^dQ 



+ 



s 

as 



JHAm 



(71) 



Here the "+" prescription with respect to (pointed in the direction of the azimuthal angle 
(h) is defined as 



2lT 



+ 00 



[/(qT)]+,qT^(qr)c?^qr = 
It acts on the asymptotic piece (69). 



[/(qr) - f{0)]g{qT)qTdqT- 



(72) 
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The functions c^^^^'^^^ are the finite residuals from the cancellation of soft singularities 
and factorization of collinear singularities. In the DREG scheme they are 



DREG: cy/'l(?) 



1 



1-0 



T 



5(1-0 



i/G 



(0 



01-0; 



2(1 -0; =u,u,d,d,... 



In the DRED scheme, these functions are 



DRED: c 



mi] 

i/j 

sum 



(0 
(0 



It can be seen from Eqs. (69) and (71) that the dependence on the the lepton decay angles 
6 and (f) in the small-gr limit is contained entirely in the overall factor G^^,^'^'^\6, 0), which 
is exactly the same as in the leading-order cross section (62). Such factorization of the 
dependence on 6 and from the hadronic dynamics refiects preservation of the chirality 
structure of the electroweak vertex as g-p ^ 0, i.e., the Born-level angular distributions of 
the leptons are not affected by the soft or collinear parton radiation. This feature, in its 
turn, automatically follows from the use of the DRED scheme in our calculation. 



(0 

(0 = 0. 



-SijCp 



71 



5(1-0; 



(73) 

(74) 
(75) 

(76) 
(77) 



B. Resummation of large logarithms 



The 0{as) cross section (71) is given in terms of generalized functions and can be used 
to predict the integrated rate. For instance, the integrated rate over the region 
(where Pt is small enough for the approximation (69) to be valid) is 



p2 



da 



[Ha,Hb) 



dq^ 



dQ'^dy dqj.dVt s 



j,k 



05/1 

TT V2 



"5 , f f^l 



p. 



(Ha) 



TT 



piHB) 
^k/b 

, AHa) 
' J a/A 



AHa) 
J a/ A 



{Hb)i 



{ H ) 

{xA,fiF)f-k/B y^B,fJ'F) + 



22 



^kb ^ Ji 



{Hb) 
b/B 



(78) 



However, Eq. (71) cannot be used yet on its own to describe the distribution over qr. To 
obtain a meaningful distribution in the small-gr region, we proceed further and identify 
Eq. (71) as an 0{as) part in the perturbative expansion of the small-gr resummed cross 
section [47], given by 



da 



(Ha,Hb) 
AB 



AB 



(g,6,y,l]) + F, 



{Ha,Hb) 
AB 



{Q,qT,y,^). 



(79) 



dQ^dydq^dn s J (27r)2' 

In Eq. (79), the first term on the right-hand side dominates when qr — ^ 0. It is expressed 
as a Fourier-Bessel transform of a form factor W^'^'^^\Q, b, y, Q) in the space of the impact 
parameter b (conjugate to qr)- 

b,y,n)= J2 E '"""V, ^) exp (-5(^-^-)(g, 6, y, Ci, C^)) 

j=u,d,... k=u,d,... 
Ci 



+ (j ^ k) 



(Hb) 
b/B 



[XB,b^, — = — 
62 0* 



(80) 



Here Ci and C2 are constants of order unity that determine the momentum scales sepa- 



rating the exponential soft factor e ^ from the convolutions 



b/A 



(x,6,Ci/C2,/i). 



The constant C3 specifies the factorization scale /i, which separates the function 



C\^l^{x,b,Ci/ C2, jji = Cs/b) from the PDF /^^^(x, = C3/6) in the convolution 



'^a/b 



fl^^ (x, 6, Ci/C2,/i = C-i/b). For definiteness, the form factor W^'^'^^\Q,b,y,Q) in 
Eq. (80) is written in its common (and, perhaps, the simplest) form proposed in Ref. [47]. 
Namely, the perturbative parts of W^'^'^'^\Q,b,y,Q) are evaluated at a variable 

(81) 



1 + byb 



2 

max 



where the parameter bmax is of order 1 GeV~\ The variable b^ is introduced 
to separate short- and long-distance dynamics in W^'^'^'^\Q,b,y,Q,). Furthermore, 
W^'^'^^\Q,b,y,^l) is normalized to the Born-level hard contribution (included in the nor- 
malization factor G^^^'^'^^\6,(f)), while higher-order corrections to the hard part are ab- 
sorbed in the functions [C®/] and . Recently discussed alternatives to this representation 
can be found in Refs. [55, 56, 89, 90]. 
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In this representation, the soft (Sudakov) factor S^^^'^'^\b,Q,y) consists of the pertur- 
bative part Sp{Q,b,Ci,C2) and the nonperturbative part Sl^p'^'^\Q,b,y): 



5(^-^-)(g, 6, y) = Sp{Q, b,, Ci, C2) + ^p"'"""^ (Q, b, y), 



where 

Sp{Q,b,Ci,C2) 



A (as (jl) , Ci) In + B (as (jl) ,^1,6-2) 



The functions A and i3 in iSp can be calculated in perturbation theory, as 



(82) 



(83) 



A\as{.^J),Cl 
B(as{]l),CuC2 



m=l 
00 



IT 



i3^-J(Ci,C2). 



m=l 



(84) 
(85) 



The coefficients in these expansions are known up to order (9(a|). Explicitly, they are 

^[^1 = Cp, 



Cp 
= 2Ci.ln 



67 TT^ 



36 12 



boC2 

B^'^ = CplCp(^-^^-3a3) 



(86) 
(87) 

(88) 



11 



36 



-TV 



193 

Is" 



+ Np{-Ln^ + ^ 
' 18 24 



/67 



vr 



iVc 

VIS 6 



-Np 
9 



In 



;C(3) 

C2bo 



+ /5o 



2 / ^0 \ . 9 ^ 3 



InMgj -ln^C2-|lnC2 



(89) 



In these equations, Np is the number of active quark flavors, [3q = (llA'^c — 2Np) /6, ({3) ~ 
1.202 is the Riemann zeta function, and 60 = 2e~'^'^ is a function of the Euler constant •yp ~ 
0.577. In the present work, we have re-derived the coefficients A^^\ i^W, and they agree with 
the earlier published results [19, 30]. The coefficient ^'^1 for the spin-dependent collisions was 
discussed earlier in [19]; since it comes from the leading-logarithmic contributions generated 
by the soft gluons, it is expected to be independent of spin.^ 



^ The spin independence of soft contributions is evident from the Feynman rules for the soft eikonal ap- 
proximation; see, for instance, Ref. [45]. 
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We have also determined the coefficient i?'^' in Eq. (89) by employing the method of 
Refs. [89-91] and order-by-order independence of the full form factor W^'^'^'^^ from the 
choice of the scales Ci/b and C2Q- By utilizing the independence on Ci and C2, we find 



(90) 



^[^1 = Si^l + /3o (4^1 In^ + In - In^ C^) - In 

where the subscript "c" denotes the coefficients in Eqs. (86)- (89) evaluated for the "canonical" 
variables Ci = bo and C2 = 1: 

411 = Cf, BW = -|Cf, (91) 
and 

-^Kt-I- ^«^>) (s^^ - f 4^ p') -^^^^ i-Ts-' S) 

To derive the last expression, we use the generalization of the relationship from [89, 90], i.e.. 



sm^.iA^ + ^„l^^ + ^y (94) 

This equation relates the canonical value of the coefficient B^'^^ to the finite contribution 
Vqq/4: = Cf{t^'^ — 8)/4 from the 0{as) virtual correction <l>(a5,£:) and coefficients 7^^ 
of the 5(1 — x) terms in the two-loop qq splitting functions on the sides of the hadrons A 
and B, respectively. Due to the helicity conservation in the quark-gluon vertex for massless 
quarks, the 5(1 — x) terms are the same in the unpolarized and polarized case, so that 

7? ^ 7|' -Cl{i- ^6C(3)) . C.N. (II . - 3C(3)) - \c.N. (i . 2^) , (95) 

Hence, similarly to J^^\ the coefficient B^'^^ in Eq. (83) is also independent of the spin of the 
incoming quarks. 

Alternatively, the spin independence of the coefficient S'^l can be shown in the anticom- 
muting 75 or the DRED schemes by analyzing the spin structure of the C(a|) cut diagrams 
with the radiation of one or two gluons. The coefficient S^^l is generated by the diagrams 
that reduce to the Born level diagram in the limit when the unobserved radiated partons 
are soft. All such cut diagrams have only one closed fermion line and can be reduced to 
the corresponding unpolarized diagram (up to an overall sign) by commutation of the 75 
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matrices from the fermion spin projectors along the fermion hne. Therefore, they lead to 
the same expression for B^"^^ as in the unpolarized case.'^ To summarize, our study confirms 
the spin independence of all Sudakov coefficients in Eqs. (86)- (89), which agree with their 
values known from the unpolarized calculations [47, 48, 92, 93]. 

Similarly to the perturbative Sudakov factor, the functions C^/^ (C, b, fi) can be calculated 
order by order in PQCD: 



m=0 ^ 



H=U,P. (96) 



The coefficients C^^^^"^^ are known up to order 0{as)- The lowest-order coefficients are 
trivial: 

47°](e,6/i) = 47°i(e,&/i) = w-O; (97) 

ClT' = ^oT = ^oT = 0; J = n, n, J, . . . . (98) 

All non-zero 0{as) coefficients in the DRED and DREG schemes are given by Eqs. (53), (54), 
and (59)-(61). As discussed in Section III, the coefficients Cj^^^^\c,,bfi) depend on the 
prescription for continuation of the 75 matrix to n 7^ 4 dimensions, and this dependence 

(H) £ 

is entirely determined by the form of the n-dimensional splitting functions P^/f, (0 the 
DRED and DREG schemes given in Eqs. (52) and (55)- (58). The above equations explicitly 
demonstrate the factorization of the collinear contributions associated with the hadron A 
from the collinear contributions associated with the hadron B. That is, the C-function 
associated with the hadron A depends only on the polarization of A and does not depend 
on the polarization of the hadron B; and vice versa. Such factorization completely agrees 
with the general structure of the resummed cross section, in which the part associated with 
the collinear radiation along the beam A is independent from the collinear radiation along 
the beam B. 

Our results for the functions C|^(^,6, /i) and C|7(^,6, /i) differ from the corresponding 
expressions in the earlier publications [19] and [30]. Firstly, in accordance with the factor- 
ization of collinear contributions, our expression for the unpolarized function cl'^^{^,b, fi) 



^ We thank W. Vogelsang for pointing out the alternative proof of the spin independence of the coeffi- 
cient ^[21. 
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is the same in the unpolarized and single-spin processes. On the other hand, Ref. [30] has 
found CflQ^\^,h, jj) in single-spin W boson production to be 

^ . = 7-n^S(ainf^), (99) 

so that the above factorization does not hold.^ This difference may be caused by the choice 
of the normalization of the q^-^^G^^^ cross section in Ref. [30]. Namely, the same result as 
in Eq. (99) would be obtained if the q^^^G^^^ cross section in n = 4 — 2e dimensions were 
averaged over two polarizations of the initial gluon and not over 2(1 — e) polarizations, as 



it is required for an unpolarized initial-state gluon in the MS scheme. 

Secondly, Refs. [19, 30] have used a different factorization scheme for the polarized 
gluon-initiated functions C^.^^{^,b, jj). In that scheme, C^^Q^\^,bj2) is identically zero for 
the "canonical" scale fi = bo/b: 

c™K.M„.,. . = -Pl2iiM^l (100) 



Refs. [19, 30] \k 



so that C^.^^^^\c,, bfi) = if /i = bo/b. To obtain such C— function, the 0{as) quark distribu- 
tion in the gluon has to be 

^^/^K-^ief. [19. 301 ^ S i-y^oii) + 1 - «) + O (4) . (101) 

which is different from /g^G(^,/i) = -(a5/27re)Pj/^(0 + C(a|) in the WS scheme. While 
this definition is not contradictory, it is not the same as the convention of the MS scheme 
[86-88]. Consequently, the functions Cf/^f^'(^, 6/i) and cf/^j'^^(^, 6/i) found in Refs. [19, 30] 
cannot be combined with the MS parton distributions. Instead, a calculation in the MS 
scheme must use the functions Cl^^^^\c,,bfi) presented in Eqs. (59)-(61). 

The remaining part of the form-factor W^'^'^^\Q , b, y, Q) to be discussed is the nonper- 
turbative Sudakov function SNp{b,Q) (see Eq. (82)), which determines the behavior of this 
form factor at 6 > 1 GeV~^. While at the present stage it is not possible to calculate this 
nonperturbative function from the first principles, its phenomenological parametrization can 
be found from the global analysis of VBP data. At the time of writing of this paper, the 
latest parametrization of SNp{b, Q) in the 6* approach in the unpolarized case is available in 



Note that Refs. [19, 30] expand in the series of (Q;s/27r), so that the 0{as) coefficients in those papers 
are two times larger than the 0{as/TT) coefficients presented here. 
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Ref. [58]. Finally, we comment on the finite term, g^, in Eq. (79). It is given by 

the diff'erence between the 0{as) finite-order cross section and its singular part (asymptotic 
piece) that is already included in the W term: 



da 



(Ha,Hb) 
AB 



dQ'^dydqpj.dVt 



da 



{Ha,Hb) 
AB 



Oias) 



dQ^dy dqj^dVL 



(102) 



C'(as),gT-*o 



The 0{as) perturbative cross section for non-zero qr is given in Eq. (67), while the asymp- 
totic piece is given in Eq. (69). This term is regular and, in fact, vanishes as qr 0, so 
that the distribution (79) in the small-gy region is approximated well by the W term. At 
It ^ Q, the cross section (79) reduces to the finite-order cross section (67), up to higher- 
order corrections. 



V. SUMMARY 



In this paper, we presented fully differential distributions for production and decay of 
electroweak bosons (7*, W^, Z^) in the collisions of proton beams of arbitrary longitudinal 
polarizations. One of our new results is the complete set of spin-dependent structure func- 
tions of order 0{as) for the angular distributions of leptons from the vector boson decay. 
This finite-order cross section is combined with the all-order sum of leading logarithms in 
the small transverse momentum region using the Collins-Soper-Sterman (impact parameter) 
resummation method [47]. The perturbative coefficients are presented at the one-loop level, 
with the perturbative Sudakov factor presented at the two-loop level. For the first time, we 
have explicitly demonstrated the universality of QCD factorization in the perturbative part 
of the spin-dependent resummed cross section. Namely, the soft factor is independent of 
spin, while the form of the "6-dependent PDFs" [Ca/b ® fb/A\{x,h, jj) is entirely determined 
by the type and polarization of the corresponding initial-state hadron, and it is independent 
from the type and polarization of the other initial-state hadron. We have also demonstrated 
the spin independence of the C(«|) coefficient B^"^"* based on the method of Refs. [89, 90], 
and, alternatively, on the analysis of the spin structure of C(a|) Feynman diagrams. The 
spin-dependent cross sections were presented in the dimensional regularization [62, 63] and 
dimensional reduction [59] schemes. 

The effects of soft parton radiation and vector boson decay considered here must be 
understood well in order to correctly describe production of Drell-Yan pairs and massive 
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electroweak bosons at the Relativistic Heavy Ion Collider. In the accompanying paper [41], 
we use these results to study measurements of sea quark parton distributions in weak boson 
production with longitudinally polarized proton beams. 
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APPENDIX: THE STRUCTURE FUNCTIONS FOR THE FINITE-ORDER 
CROSS SECTION 

This appendix collects explicit expressions for the structure functions ''TabiH^, Hb, 's, t, u) 
that enter the 0{as) real emission cross section in Eq. (67). Eqs. (104-110) show the 
independent structure functions '^Tgg, ^Tcq, , and, in the case of the single-spin cross section, 
'^TqG- The end of the appendix includes the rules to derive the rest of nonzero structure 
functions. As a convenient notation, the mathematical expressions will contain the functions 
T±{a,b), defined as 



The independent structure functions in the unpolarized cross section are as follows [50]. 

In the q'^^\pa)(f'^\pb) subprocess, 




(103) 



-1 




0. 



Irji 



T+{u,t) 
Ml 



J- n 



Q T_{u, t) 
qr Ml ' 
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4„ T_{u,t) 



^ _ T+(g,M) t 



In the G'^^\pa)q^^\ph) subprocess, 
Or - - _T±(zl2^i 

^Gq iGg - 

1 ^ Q T4u,?) + iQ'-uf t 
3 ^ g T+(n,s)-2n(g^-g-) f 

'^'^ Mt r 

. T+(u,s) +2(Q^ -s)st , , 

Next, the independent structure functions in the single-spin cross section are as follows. 

In the q^^\pa)q^^\pb) subprocess, 
_i _ T+{u,t) 

~ M| ' 

~ qr ' 

J- ( 



11 ~ /\/r_ ^2 



Mt qT 



In the q'^^\pa)G^^\pb) subprocess, 
-It. _ T+(t,s)u 

J-qG — 2 ~' 

qr s 

, _ Q T_{%u) + {Q^-tf u 
qr M| r 

3 ^ Q T+{t,s)-2iQ'-s)t u 
Mt ql r 

= M^T f ^ 
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In the G'^^\pa)q^^\pb)suhpvocess, 

J^Gq — 2 ~' 

ot - 2t - ^-(^' ~^) ^ 

1 _ Q T4uS)-iQ'-u)n 

3 ^ Q T_iu,s) + 2uiQ'-s) t 

""^ Mt qi r 

4 ^ T_{u,s)-2s{Q'-s) t 
= ^ M^^? f 

Finally, the independent structure functions in the double-spin cross section are as follows. 

In the q^^\pa)q^^\pb) subprocess, 
_i _ T+{u,t) 

M| ' 
1 _ Q T_(u,?) 

qr Ml ' 



3 



11 T\/r_ ^2 



Mt q^ 
4^ T- (u, t) 

In the G^^\pa)q^^\pb) subprocess, 



"^^Gq -2 



T_ (s', u) t 

2 ~ ' 



Ot^ _ 2t. _ T_{-s,u) t 
T,G - To,- 

1 ^ g T„(^,t)-(g^-72)^ ? 

"""^ gT M| s' 

3 ^ g T_(t2,g-) + 2^(Q^-g) f 
Mt ql s' 

= -2 ^ \^ 110 

All other non-zero structure functions can be obtained from Eqs. (104-110) by applying 
the following transformation rules: 
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p = - 


-1,0,1,2 


P 


= 3,4 


PT- 


= ±"^.5 


PT- 
qq 






= ±'T,G 




= T'T^G 




= ±'TGq 


"Tcq 


= T'Tcq 



The upper sign in these rules should be used to obtain the structure functions in the 

unpolarized and double-spin cross sections. The lower sign should be used to obtain the 

structure functions in the single-spin cross section. An additional rule should be applied to 

the unpolarized and double-spin cross sections in order to relate the structure functions in 

the qG and Gq subprocesses: 

P = -l,0,2,4 P = l,3 
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